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Although microbial communities have varying degrees of exposure to environmental 
stresses such as chemical pollution, little is known on how these communities 
respond to environmental disturbances and how past disturbance history affects 
these community-level responses. To comprehensively understand the effect of 
organophosphorus insecticide application on microbiota in soils with or without 
insecticide-spraying history, we investigated the microbial succession in response to 
the addition of fenitrothion [0,0-dimethyl 0-(3-methyl-p-nitrophenyl) phosphorothioate, 
abbreviated as MEP] by culture-dependent experiments and deep sequencing of 16S 
rRNA genes. Despite similar microbial composition at the initial stage, microbial 
response to MEP application was remarkably different between soils with and without 
MEP-spraying history. MEP-degrading microbes more rapidly increased in the soils with 
MEP-spraying history, suggesting that MEP-degrading bacteria might already exist at a 
certain level and could quickly respond to MEP re-treatment in the soil. Culture-dependent 
and -independent evaluations revealed that MEP-degrading Burkholderia bacteria are 
predominant in soils after MEP application, limited members of which might play a 
pivotal role in MEP-degradation in soils. Notably, deep sequencing also revealed that 
some methylotrophs dramatically increased after MEP application, strongly suggesting 
that these bacteria play a role in the consumption and removal of methanol, a harmful 
derivative from MEP-degradation, for better growth of MEP-degrading bacteria. This 
comprehensive study demonstrated the succession and adaptation processes of microbial 
communities under MEP application, which were critically affected by past experience of 
insecticide-spraying. 
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INTRODUCTION 

Biological communities are exposed to varying levels of envi- 
ronmental stresses or disturbances such as global warming, 
typhoon, drought, and bush fire (Phillips et al., 2009; O'Connell 
and Nyman, 2011; Penuelas et al, 2013; Rota et al, 2014), 
and little is known about community succession and adapta- 
tion processes under such environmental stimuli. In addition 
to naturally occurring environmental disturbances, chemical 
insecticides have been developed and used worldwide to con- 
trol agricultural, hygienic, and household pest insects. Although 
insecticides have revolutionized modern agriculture and sanita- 
tion in terms of pest management, abuse of usage sometimes 
causes serious problems including environmental pollution from 
its residues, health hazard, unexpected effects on non-targeted 
organisms, and evolution of insecticide-resistant insects (Whalon 
et al, 2008; Diez, 2010). Thus, insecticide-spraying is thought 
to be a man-made, yet strong, environmental stress on natural 
organisms. 



Organophosphorus insecticide, a general name referring to 
insecticides containing a phosphoester bond, is one of the 
most widely used insecticides. It includes a number of com- 
mercially available chemicals such as diazinon, malathion, and 
dichlorvos (Singh and Walker, 2006; Singh, 2009). The chemical 
compounds found in organophosphorus insecticide show high 
mammalian toxicity as well as insecticidal activity by inhibit- 
ing acetylcholinesterase (AchE), causing overstimulation of the 
cholinergic synapses by the overaccumulation of acetylcholine 
(da Silva et al., 2013). Fenitrothion [0,0-dimethyl 0-(3-methyl- 
p-nitrophenyl) phosphorothionate, abbreviated as MEP] is one 
of the organophosphorus compounds commonly used because 
of its broad-spectrum activity and less mammal toxicity. In 
natural fields, MEP can be degraded photochemically, but is 
mainly degraded by soil bacteria. Biodegradation pathways of 
MEP were first investigated in forest soils using C 14 -labeled 
MEP (Spillner et al., 1979) and are currently described in the 
International Programme on Chemical Safety IPCS (1992). In 
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microorganisms, MEP is hydrolyzed to 3-methyl-4-nitrophenol 
and then finally decomposed to CO2 through a phenol ring cleav- 
age pathway (IPCS, 1992). Previous culture-dependent studies 
reported that members of the genera Burkholderia, Pseudomonas, 
Cupriavidus, Corynebacterium, Arthrobacter, and Sphingomonas 
have been identified as MEP-degraders, some of which use the 
compound as a sole carbon source (Tago et al, 2006; Zhang et al., 
2006; Kim et al., 2009). Considering that only a tiny fraction 
of environmental bacteria is culturable (Wilson and Piel, 2013), 
culture-independent as well as culture-dependent approaches are 
pivotal to comprehensively understand microbial dynamics in 
soils after spraying with MEP, as previously mentioned (Jacobsen 
and Hjelms0, 2014). 

Recent advances in sequence technologies, the so-called 
next generation sequencing (NGS), provide a faster and sim- 
pler alternative to previously established techniques for the 
comprehensive investigation of microbial community struc- 
tures (Moorthie et al., 2011). Fairly recent works using NGS 
include the evaluation of microbial community compositions 
of soils polluted with organochlorine pesticides including 
Dichlorodiphenyltrichloroethane (DDT), hexachlorocyclohex- 
ane (HCH), and 2-chloro-4-ethylamino-6-isopropylamino-s- 
triazine (ATZ) (Fang et al, 2014). Based on the metagenomic 
survey of functional genes, these studies have suggested that 
diverse bacterial species complement one another to degrade 
these organochlorides (Fang et al., 2014). 

In crop fields, biodegradation of the insecticides is one 
of the pivotal factors affecting the recovery and sustainabil- 
ity of soil ecosystems. Furthermore, it has been reported 
that one or more previous applications of the same insecti- 
cides and/or other structurally-analogous insecticides frequently 
result to less persistence of the insecticide in the environment 
(Felsot, 1989). In soil environments where the same insecti- 
cide has been continuously sprayed, it has been hypothesized 
that insecticide-degrading microbes grow by assimilating the 
insecticide and are enriched in the soil, which cause exces- 
sively enhanced biodegradation of the sprayed insecticides 
(Arbeli and Fuentes, 2007). However, it remains unclear how 
microbial communities transit under and adapt to insecticide- 
spraying and how the history of insecticide-spraying affects 
these succession and adaptation processes of soil microbial 
communities. 

By combining the culture-independent approach using 
the NGS technique with culture-dependent isolation of 
MEP-degrading bacteria, we investigated the detailed process 
of microbial succession when MEP was sprayed on field- 
collected soils with or without a history of MEP-spraying 
prior to collection. Among the collected soils, one has been 
sprayed with MEP at least for 4 years and the other has never 
been sprayed with MEP nor with other organophosphorus 
compounds. The NGS analysis combined with culture- 
dependent approach clearly demonstrates previously unseen 
dynamics of soil microbiota during MEP-spraying, high- 
lighting the different responses of microbiota against the 
experimental MEP-spraying in the different types of soil, 
and identifying the key players of MEP-degradation in soil 
environments. 



MATERIALS AND METHODS 
SOILS 

We collected two andisol soils from agricultural fields located 
at the National Institute for Agro-Environmental Sciences 
(Tsukuba, Ibaraki, Japan; 36°0'N, 140° l'E): soil S (Sprayed soil) 
was collected from a crop field that had been sprayed with MEP 
for at least 4 years [pH(H 2 0), 6.3; water content (w/w), 35.2%; 
total carbon, 5.0%; total nitrogen, <0.3%] ; soilN (Naive soil) was 
collected from a field that had never been sprayed with MEP and 
other organophosphorus compounds [pH(H20), 6.6; water con- 
tent (w/w), 33.7%; total carbon, 4.4%; total nitrogen, <0.3%]. 
Characterization of these soil samples is summarized in Table SI. 
Each soil was passed through a 2 mm sieve to remove large organic 
litters like leaves and roots. 

INSECTICIDE-SPRAYING EXPERIMENTS 

Experimental design of the insecticide-spraying test is shown in 
Figure SI. Approximately 150 g (dry weight) of each of the sieved 
soils were transferred into plastic pots (10 x 8 cm; opening diam- 
eter x depth). The potted soils were incubated at 25° C with 
weekly application of MEP solution (diluted by distilled water) 
at a specific loading of 50 mg kg _1 -dry soil (Tago et al., 2006). 
For both soils S and N, duplicate runs (pot 1 and pot 2) were 
prepared and individually subjected to experiments and further 
analyses. Control treatments were conducted under the same con- 
ditions, but the addition of MEP solution was replaced by equal 
amount of distilled water only once a week. These pots were cov- 
ered with a sheet of aluminum foil during incubation to prevent 
water evaporation and photodegradation of MEP. Soil samples 
from the surface layer up to a depth of 1 cm were collected 1 week 
after the final spraying, e.g., soil samples treated three times with 
MEP were collected 1 week after the 3rd MEP-spraying (see Figure 
SI), and then subjected to the various analyses described below. 

CFU COUNTING OF MEP-DEGRADING BACTERIA 

The total CFU (colony forming units) counts of MEP-degrading 
bacteria was measured by plating serial dilutions of soil sam- 
ples [1% (w/w) diluted in distilled water] on 1.5% agar plates 
of MEP medium [0.08% MEP, 20 mM potassium phosphate (pH 
7.0), 0.1% (NH 4 ) 2 S0 4 , 0.02% NaCl, 0.01% MgS0 4 .7H 2 0, 0.05% 
CaCl 2 .2H 2 0, 0.0002% FeS0 4 .7H 2 0, 0.1% yeast extract] and 
incubated for 4 days at 25° C, as previously described (Tago et al., 
2006). The CFU of MEP-degrading bacteria were determined by 
counting the colonies with a halo on the MEP plates. In the halo- 
forming assay, the bacteria that are capable of degrading MEP to 
3-methyl-4-nitrophenol were detected. 

ISOLATION AND IDENTIFICATION OF MEP-DEGRADING BACTERIA 

Isolated colonies were identified based on the partial sequences 
(c.a. 600 bp) of their 16S rRNA genes. During isolation, each 
of the colonies was picked with a sterile toothpick and then 
suspended in 35 [i\ TE buffer [10 mM Tris-HCl, 1 mM EDTA 
(pH 8.0)] in each well of a 96 well titer plate. The cells 
were subjected to heat shock treatment (98°C for 3min) in 
order to extract the genomic DNA, which was used as tem- 
plate for PCR amplification. The PCR amplification of bacte- 
rial 16S rRNA gene was performed with the use of a KOD Fx 
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Neo polymerase (TOYOBO, Tokyo, Japan), a universal primer 
set, 16SA1 [ 5' - AGAGTTTGATCMTGGCTCAG- 3' ] and 16SB1 
[ 5' -TACGGYTACCTTGTTACGACTT- 3' ] (Fukatsu and Nikoh, 
1998), and the template DNA extracted as described above. 
The temperature program for PCR is as follows: 94° C for 
2min, followed by 20 cycles of 98°C for 10 s, 55°C for 30 s, 
and 68°C for 90s. The resulting 1.5 kb amplicons were used 
as templates of sequencing reaction conducted with a uni- 
versal primer 357F [5'-CCTACGGGAGGCAGCAG-3'] (Muyzer 
et al, 1993). Taxonomic assignment of the resulting sequences 
was double checked by the RDP multiclassifier ver. 1.1 (Wang 
et al., 2007) and BLASTN-search (http://ncbi.nlm.nih.gov/blast/) 
against GreenGene database (DeSantis et al., 2006). 

DEEP SEQUENCING 

DNA was extracted from soil S and soil N samples before 
treatment (SO, NO), after the 2nd (S2) and 3rd (S3, N3) 
MEP-treatments, and after the 3rd control treatment (S3C, 
N3C) by using a slightly modified published protocol (Noll 
et al, 2005). Briefly, 0.5 g wet soil was combined with 0.08 g 
skim milk, 700 u.1 extraction buffer [50 mM Tris-HCl (pH 
8.0), 1.7% Polyvinylpyrrolidone, 20 mM MgCl 2 ], and 0.5 ml 
zirconia/silica beads (0.1mm diameter, Biospec, OK, USA), 
and blended in a Multi-Beads Shocker (Yasui Kikai, Osaka, 
Japan) at 2000 rpm for 60s. Crude extracts were purified 
by phenol/chloroform/isoamyl alcohol (25:24:1) treatment fol- 
lowed by isopropanol precipitation. Contaminating RNA was 
digested with Ribonuclease (Nippon Gene, Toyama, Japan). 
The prepared DNA was subjected to PCR amplification of 
16S rRNA gene for deep sequencing. The variable region (V4) 
of bacterial 16S rRNA gene was amplified using universal 
primers 515F [5'-GTGCCAGCMGCCGCGGTAA-3'] and 806R 
[5'-GGACTACHVGGGTWTCTAAT-3'] (Caporaso et al, 2012). 
The PCR reaction mixture comprised of 50 |xM each dNTP, 
0.4 |xM 515F with Illumina P5 sequences, 0.4 |xM 806R with 6- 
bases indexes and Illumina P7 sequences (Caporaso et al., 2012) 
(Illumina, San Diego, CA, USA), Q5 High-Fidelity DNA poly- 
merase with Q5 reaction buffer (New England BioLabs, Ipswich, 
MA, USA) and the extracted soil DNA template. The PCR con- 
ditions were as follows: initial denaturation at 98° C for 90s, 
followed by 25 cycles of 98°C for 10s, 54°C for 30s, 72°C for 
30 s, and a final extension at 72°C for 2min. The PCR ampli- 
cons were purified as described previously (Itoh et al., 2014). 
DNA libraries containing all tagged-amplicons and internal con- 
trol phiX were generated for paired-end sequencing by the MiSeq 
sequencer using MiSeq Reagent kit v2 (Illumina, San Diego, CA, 
USA) according to the manufacturer's instruction. 

QUANTITATIVE PCR 

Quantitative PCR (qPCR) of bacterial 16S rRNA gene was per- 
formed to amplify bacterial 16S rRNA genes using a Power SYBR 
Green PCR Master Mix (Applied Biosystems, Foster City, CA, 
USA) and the LightCycler 96 System (Roche Applied Science, 
Indianapolis, IN, USA). The reaction mixture was comprised 
of 2x SYBR Green PCR Master Mix, 0.2 [iM 515F and 806R 
primer pairs (Caporaso et al., 2012), 0.5(ig/[il BSA, and soil 
DNA as a template. The PCR conditions were as follows: initial 



denaturation at 95°C for lOmin, followed by 45 cycles of 95°C 
for 30 s, 57° C for 30 s, and 72°C for 30 s. The amount of bac- 
terial 16S rRNA gene copies was calculated on the basis of the 
standard curve constructed using a dilution series of the tar- 
get PCR product of Burkholderia sp. SFA1 (DDBJ accession no. 
AB232333). 

DATA ANALYSIS 

Internal control phiX sequences and low-quality sequences were 
removed, and paired sequences were joined as described previ- 
ously (Itoh et al., 2014). Chimeric sequences were removed using 
uchime algorithm on Mothur program ver. 1.29.2 (Schloss et al., 
2009) through the alignment with Greengene database (DeSantis 
et al, 2006). The resulting sequences were subjected to taxonomic 
assignment by using RDP multiclassifier ver. 1.1 (Wang et al., 
2007) with an 80% confidence threshold. Operational taxonomic 
units (OTU) based analyses, including estimation of diversity 
indices and Principal Coordinate Analysis (PCoA) based on OTUs 
with 3% differences, were performed by using the macqiime ver. 
1.6.0 (Caporaso et al., 2010). Phylogenetic tree was constructed by 
the neighbor-joining method with the bootstrap test (1000 repli- 
cates) under MEGA ver. 4.0.2 (Tamura et al., 2007). Analysis of 
variance (ANOVA) was performed using R software ver. 3.0.1 (R 
Development Core Team, 2008) to analyze differences in qPCR 
data among samples. 

NUCLEOTIDE SEQUENCE ACCESSION NUMBER 

The nucleotide sequences reported in this study were deposited 
in the DDBJ/Genbank/EBI databases under the accession num- 
bers: AB904935-AB905196 (16S rRNA gene sequences of isolates, 
Table S2) and the MG-RAST database (http://metagenomics.anl. 
gov/, Meyer et al., 2008) as a "Microbial succession under MEP- 
spraying in 2012" project under the ID: 4562358.3-4562369.3 
(deep sequencing, Table 1). 

RESULTS 

POPULATION DYNAMICS OF MEP-DEGRADING BACTERIA IN 
MEP-SPRAYED SOILS 

Abundance changes of MEP-degrading strains by a culture- 
dependent method were investigated using two types of soils 
with different MEP-spraying histories: MEP-sprayed "soil S" and 
MEP-naive "soil N" (Table SI). Although the CFU of MEP- 
degrading bacteria were below the detectable level (<10 3 cfu 
g _1 dry soil) even 1 week after the insecticide application, 
they became detectable after the 2nd and 3rd treatments in 
soil S and soil N, respectively, (Figure 1). In the control, with 
only the addition of distilled water, no MEP-degraders were 
detected even after the 3rd treatment (Figure 1). These results 
indicate that the repeated spraying of MEP acts as a strong 
selective pressure on soil microbiota, confirming our findings 
from previous studies (Tago et al., 2006; Kikuchi et al, 2012). 
It has also been reported that some bacteria are capable of 
utilizing MEP as a sole carbon source, fostering their popu- 
lation increase (Hayatsu et al, 2000; Tago et al, 2006; Zhang 
et al., 2006; Kim et al., 2009), which likely explains the above- 
observed population dynamics of MEP-degrading strains in soil 
environments. 
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Table 1 | Summary of deep sequencing and qPCR. 



Soil 


No. of 
spraying 


Pot No. 


Libarary 
ID 


No. of 
sequences 6 


No. of 
OTU5.M 


Cx d 


Diversity indices' 9 
chaol shannon 


1/simpson 


No. of 

copies' 




S 


0 


1 


S0p1 


13,202 


1,206 


0.97 


1,944 ± 107 


7.98 ±0.04 


65±2.7 


3.6 ±0.1 x 


10 8 


s 


2 


1 


S2p1 


13,346 


796 


0.99 


1,188 ± 103 


5.9 ±0.03 


15 ±0.3 


4.5 ±0.3 x 


10 s 


s 


3 


1 


S3p1 


14,600 


575 


0.99 


813 ± 79 


3.86 ±0.05 


4.5±0.1 


4.7 ±0.2 x 


10 s 


s 


0 


2 


S0p2 


6,218 


1,062 


0.93 


2,446 ± 78 


8.62 ±0.01 


109 ±0.7 


4.6 ± 0.2 x 


10 s 


s 


2 


2 


S2p2 


10,223 


1,154 


0.94 


2,324 ± 116 


7.72 ±0.03 


42 ±0.8 


4.0 ± 0.6 x 


10 s 


s 


3 


2 


S3p2 


8,213 


897 


0.96 


1,731 ± 100 


6.43 ±0.03 


13±0.3 


4.4 ± 0.5 x 


10 s 


s 


0 a 


control 


S3C 


17,933 


1,533 


0.94 


2,260 ± 123 


8.6 ±0.04 


124±2.8 


4.8 ± 1 .2 x 


10 s 


N 


0 


1 


N0p1 


13,540 


1,500 


0.85 


3,268 ± 190 


9.07 ±0.04 


138±5.8 


5.5 ± 1.8 x 


10 s 


N 


3 


1 


N3p1 


10,171 


1,028 


0.95 


1,919 ± 114 


7.3 ±0.03 


36 ±0.8 


6.3 ± 0.9 x 


10 s 


N 


0 


2 


N0p2 


15,054 


1,747 


0.87 


3,655 ± 180 


9.1 ±0.02 


151 ±3.4 


6.3 ± 0.2 x 


10 s 


N 


3 


2 


N3p2 


34,865 


1,970 


0.90 


2,977 ± 187 


7.78 ±0.04 


33 ±1.2 


6.8 ± 0.5 x 


10 s 


N 


0 a 


control 


N3C 


72,066 


2,459 


0.90 


3,470 ± 139 


9.06 ±0.04 


149±9.3 


8.6 ± 1.0 x 


10 8 



a Control pots received distilled water (DW) once a week: soils were collected 1 week after 3rd DW-spraying. 
b Number of sequences after removal of low-quality, chimeric, and archaeal sequences. 
0 Clustered at the 3% distance level. 

d Calculated from the equation, Cx= 1-(n/N), where "n" is the number ofOTUs composed of only one sequence (singleton) and N is the total number of sequences. 
"Each index was calculated based on the same amount of sequences (5174) sub-sampled from original libraries (Mean ± SDs of 10 time sub-samplings is shown). 
1 Copy numbers (copies g~ ] -dry soil) of 16S rRNA genes estimated by qPCR (Mean ± SDs of three time measurements of the same DNA sample is shown). 
Analysis of variance (ANOVA) was performed using R software ver. 3.0. 1 (R Development Core Team, 2008), showing the copy numbers were not significantly 
different between samples that were originated from the same field soil. 
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FIGURE 1 | Transition of CFU of MEP-degrading bacteria in MEP-sprayed soils. CFU of MEP-degrading bacteria in soil S (A) and soil N (B) were counted by 
using MEP medium plates. Results of replicated experiments (pot 1 and pot 2) are shown. ND, not detected (< 1 .0 x 10 3 cfu g _1 -dry soil). 



ISOLATION AND IDENTIFICATION OF MEP-DEGRADING BACTERIA 

From the soil S after the 2nd and 3rd MEP- treatments (S2 and 
S3), and soil N after the 3rd MEP-treatment (N3), colonies of 
the MEP-degrading bacteria were isolated and identified based on 
the partial sequences of 16S rRNA gene. The sequence analyses 
demonstrated that the majority belonged to the betaproteobac- 
terial genus Burkholderia. Specifically, they occupied around 
97.9, 96.9, and 61.1% of the total sequenced-bacteria isolated 
from the S2, S3, and N3 soils, respectively, (Figure 2; Table S2). 
Notably, in contrast to the predominance of Burkholderia strains 
in the S2 and S3 samples, MEP-degraders from the N3 sample 
were more phylogenetically diverse with the presence of strains 
from the genera Dyella, Ralstonia, Pandoraea, Achromobacter, and 
Cupriavidus. 



SUCCESSION OF MICR0BI0TA IN MEP-SPRAYED SOILS AS REVEALED 
BY DEEP SEQUENCING 

For comprehensive understanding of the temporal dynamics 
of microbiota in the MEP-sprayed soils, we performed deep 
sequencing of bacterial 16S rRNA genes from both soils before 
treatment (SO and NO), after the 2nd (S2) and 3rd (S3 and N3) 
MEP-treatments, and after the 3rd control-treatment (S3C and 
N3C). The sequencing of PCR amplicons of partial 16S rRNA 
genes produced > 2 x 10 5 sequences in total (Table 1). The qPCR 
data showed that the amount of bacterial 16S rRNA genes was sta- 
ble at around 4 x 10 8 and 6 x 10 s g~'-dry soil in soil S and soil 
N, respectively, during the treatments (Table 1). Meanwhile, in 
both soils S and N, diversity indices, chaol, Shannon, and recip- 
rocal simpson decreased after MEP-spraying (Table 1), suggesting 
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unclassified 

Cupriavidus 
■ Ralstonia 

Achromobacter 

Dyella 
m Pandoraea 

Burkholderia 



(94) (96) (72) 

FIGURE 2 | Genus-level compositions of MEP-degrading bacteria 
isolated from MEP-sprayed soils. MEP-degrading isolates from soil S after 
two times and three times MEP spraying (S2 and S3, respectively) and 
from soil N after three times MEP spraying (N3) were identified by partial 
(c.a. 600 bp) sequencing of 1 6S rRNA gene. Numbers in parentheses 
indicate total numbers of inspected MEP-degrading strains that were 
randomly selected from 2 pots in each soil (i.e., results of pots 1 and 2 are 
merged here). Sequences were subjected to taxonomic assignment by 
using RDP multiclassifier with a threshold level of 80%. 



that MEP treatment caused the decrease of soil microbial 
diversity. 

Similarities among the sequence libraries of each soil sam- 
ple were analyzed by PCoA based on weighted unifrac distance 
(Figure S2). The results suggest that the bacterial compositions 
of S2, S3, and N3 were distinct from those of SO and NO sam- 
ples, although the shift direction was different between the soil 
S and soil N. The controls, S3C and N3C, were closely related 
to the soils before treatments, SO and NO. These results also 
clearly demonstrate that MEP-spraying dramatically altered the 
community structure of soil microbiota. 

BACTERIAL GROUPS RESPONDING TO MEP-SPRAYING 

Figures 3A,B show the phylum-level (class-level in the 
Proteobacteria) distribution of sequences obtained from the 
deep sequence libraries. Among the diverse bacterial groups, 
Betaproteobacteria drastically increased after MEP-spraying in all 
samples from both soils S and N. Especially in soil S, the relative 
abundance of Betaproteobacteria was over 40% (80.4% in pot 1 
and 49.8% in pot 2) after the 3rd MEP-treatment (Figure 3A). In 
the control samples (S3C and N3C), the bacterial composition 
appeared to be stable based from the initial states (SO and NO), as 
also supported by the PCoA profile (Figure S2). 

Within the increasing community of Betaproteobacteria, the 
genus Burkholderia was the most abundant, followed by Ralstonia 
and Pandoraea (Figures 3C,D). Before MEP treatment, few 
Burkholderia sequences were detected. However, after the 3rd 
treatment, the frequencies of Burkholderia drastically increased 
up to 53.8% in soil S pot 1, 41.1% in soil S pot 2, 55.7% in soil 
N pot 1, and 66.9% in soil N pot 2. 



Figure 4 shows increase in the relative proportion of spe- 
cific bacteria genera during the MEP-spraying experiments. It 
depicts the top 10 genera with the higher frequency after the 
3rd treatment. From the list, the Burkholderia, Cupriavidus, 
Pseudomonas, and Arthrobacter are already well-known MEP- 
degrading strains (Tago et al, 2006; Zhang et al, 2006; Kim 
et al., 2009). Furthermore, the genera Pandoraea, Ralstonia, 
and Dyella, which were additionally identified earlier as MEP- 
degrading group (Figure 2), were also included in the top 10 
genera that highly responded to MEP (Figure 4). In both soils S 
and N, the frequencies of Burkholderia and Ralstonia drastically 
increased in response to MEP-spraying at two to three orders of 
magnitude, suggesting that these genera might possess a strong 
MEP-assimilation ability to highly adapt to the biological niche 
resulting from MEP-sprayed soils. Other bacterial groups, such as 
Rhodanobacter, Nevskia, and Methylobacillus, were also included 
in the top 10 (Figure 4), although they were not detected through 
the culture-dependent approach (Figure 2). These strains might 
assimilate the intermediates of MEP degradation such as 3-methl- 
4-nitrophenol and p-nitrophenol (Hayatsu et al., 2000; Arora 
et al., 2013), and thus proliferate under the MEP treatment. The 
frequencies of these bacterial groups with high increasing ratios 
after MEP treatment were almost stable in the control set-ups 
in both of soils S and N (Fold: 0.2~3.4), emphasizing that these 
bacteria respond to MEP-spraying. 

BURKHOLDERIA STRAINS RESPONDING TO MEP-SPRAYING 

To clarify the diversity of Burkholderia strains in the MEP-sprayed 
soils, 16S rRNA gene sequences of Burkholderia derived from 
the MEP-degrading isolates and the deep sequencing libraries of 
the soils treated with MEP three times (i.e., S3 and N3 soils) 
were classified into OTUs defined by 100% sequence identity. 
In the MEP-degrading isolates, 4 and 8 OTUs were identi- 
fied from the S3 and N3 soils, respectively, in which a single 
and distinct OTU dominantly occupied at > 70-90% frequency 
(Figures 5A,B). Although > 100 OTUs were determined from all 
of the S3pl, S3p2, N3pl, and N3p2 libraries, a single and distinct 
OTU was dominant with >40% frequency in each soil sample 
(Figures 5C-F). In sequence libraries of both soil S and soil N, 
the dominant OTUs were not detected or a few sequences were 
detected before the MEP-treatment (Figure S3). These culture- 
dependent and -independent studies strikingly demonstrated that 
a particular phylotype of Burkholderia dominantly responds in 
soil environments in the presence of MEP. 

We then estimated the phylogenetic relationship of these dom- 
inant Burkholderia strains based on 255 bp sequences of the V4 
region of 16S rRNA gene (Figure 6). Regardless of whether the 
strain was identified by the culture-dependent or -independent 
method, each of the dominant Burkholderia OTUs derived from 
isolates and deep sequences was clustered together with the OTUs 
detected from the soil of the same origin (Figure 6), strongly 
indicating that our isolates are representative strains of MEP- 
degrading Burkholderia dominating the MEP-sprayed soils. This 
also indicates that the Burkholderia frequently detected from 
deep sequencing in MEP-sprayed soils are most likely to have 
MEP-degrading activity. The most dominant OTUs in soil S sam- 
ples were related to pathogenic Burkholderia such as B. cepacia, 
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FIGURE 3 | Structural changes of microbial communities in 
MEP-sprayed soils. Soil S (A,C) and soil N (B,D) were analyzed by 
deep sequencing of 16S rRNA gene. For each soil, results of replicated 
experiments (pot 1 and pot 2) are shown. Relative distributions of the 
sequences at (A,B) phylum-level (class-level in the Proteobacteria) and 



(CD) genus-level within Betaproteobacteria are shown. In (CD), total 
numbers of betaproteobacterial sequences excluding those of 
unclassified genus are depicted on the bars. Sequences were subjected 
to taxonomic assignment by using RDP multiclassifier with a threshold 
level of 80%. 



B. glumae, and B. pseudomallei of the Burkholderia cepacia com- 
plex (Bcc) group (Coenye et al., 2001), while those in soil N sam- 
ples were related to known MEP-degrading bacteria isolated from 
soils and symbiotic Burkholderia in stinkbugs (Tago et al, 2006; 
Kikuchi et al, 201 1). These results suggest that the phylogenic lin- 
eage of Burkholderia responding to MEP-spraying appears likely 
to not be general, but specific for each soil environment. 

DISCUSSION 

EFFECTS OF FIELD-USE HISTORY ON SOIL MICROBIOTA RESPONDING 
TO MEP TREATMENT 

Deep sequencing combined with culture-dependent approach 
clearly demonstrated that soil microbiota is strongly affected by 
the application of MEP, and furthermore, bacterial communities 



responding to MEP-spraying were remarkably different between 
the two andisol soils (soil S and soil N) (Figures 3, 4, 6 and Figure 
S2) despite their almost similar origin and chemical properties 
(Table SI). Notably, the two soils have different histories of insec- 
ticide application: the agricultural field where soil S was obtained 
had experienced at least 4 years of MEP spraying; the field where 
soil N was collected was not subjected to such insecticide applica- 
tion (Table SI). Therefore, in soil S, MEP-degrading bacteria (i.e., 
Burkholderia species) might already exist at a certain level, thus 
allowing them to respond quickly to the presence of the chem- 
ical. Meanwhile, the greater variety of MEP-degrading bacteria 
in soil N may be due to the absence of selection pressure from 
MEP. Although the enhanced biodegradation caused by repeated 
treatment with the same insecticide were reported from the 1970s 
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FIGURE 4 | Temporal changes of relative proportions of bacterial 
genera during MEP-spraying. (A,B) Soil S; (CD) soil N. Relative 
proportions of bacterial genera were estimated based on the read 
numbers of deep sequencing of 16S rRNA sequences. Closed asterisks 
indicate the genera, which have been reported to degrade MEP (Tago 



et al., 2006; Zhang et al., 2006; Kim et al., 2009), or are newly 
identified as MEP-degrading bacteria in this study (Figure 2). Open 
asterisks indicate methylotrophs. In addition to the relative increasing 
rates in MEP-sprayed soils, those in DW-sprayed soils (controls) are 
shown in text. 



(Felsot, 1989), little is known about the relation between soil 
microbial succession/adaptation and insecticide-spraying history. 
This is the first report to show the different microbial-community 
responses against insecticide re-treatment in soils that previously 
experienced such chemical application vs. naive soil with no such 
treatment prior to the experiment. 

It should be noted that the bacterial composition of the ini- 
tial states in both soils S and N was similar to each other 
(Figure S2) despite their different MEP-spraying history. This 
result strongly suggests that, although speculative, microbial 
community-structure disturbed by insecticide-spraying could 
recover to an initial community composition once the insecti- 
cide application is terminated. To clarify the plasticity of micro- 
bial community in more detail, it would be of great interest 
to investigate succession process of insecticide-enriched micro- 
biota after terminating insecticide spraying. Interestingly, our 
results demonstrated that soil diagnostics of biological activity 
should not be inferred only based on a snapshot of environmental 



microbial community. To estimate potentials of a soil (micro- 
biota) more precisely, it may be more important to trace the 
microbial succession under environmental stimuli. 

BURKHOLDERIA IS A KEY PLAYER IN MEP-DEGRADATION 

Among the bacterial groups that dominantly responded to this 
insecticide, the genus Burkholderia was the most abundant after 
MEP-treatment (Figures 3C,D) and included a number of MEP- 
degrading isolates (Figure 2). Previous studies from China, Japan, 
and South Korea that were based on culture-dependent meth- 
ods have repeatedly isolated MEP-degrading Burkholderia strains 
from MEP-contaminated soils at a higher frequency relative to 
other genera (Tago et al., 2006; Zhang et al, 2006; Kim et al., 
2009). Taken together, our findings from both culture-dependent 
and -independent approaches strongly suggest that Burkholderia 
is a key bacterial group for MEP-degradation in soil environ- 
ments. Remarkably, Burkholderia species were also found to 
degrade other organophosphorus insecticides, which includes 
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FIGURE 5 | Relative abundance of Burkholderia OTUs in MEP-sprayed 
soils. Relative abundance of 16S rRNA gene sequences assigned to 
different OTUs (100% sequence identity thresholds) of the genus 
Burkholderia in the MEP-degrading isolates (A,B) and the deep sequence 



libraries (C-F). (A,C,E) Soil S; (B,D,F) soil N. Data showed the top 4 or 10 

OTUs with high frequency in Burkholderia isolates or sequences, 
respectively. OTUs with asterisks were subjected to phylogenetic analysis 
(see Figure 6). 



parathion, methyl parathion, glyphosate, and chlorpyrifos- 
methyl (Keprasertsupa et al., 2001; Kuldinsky-Sobral et al., 2005; 
Kim et al., 2007b, 2009). Considering that MEP-degrading strains 
of Burkholderia frequently show cross-acclimation against other 
organophosphorus compounds (Hayatsu et al., 2000; Kim and 
Ahn, 2009; Kim et al, 2009; Kikuchi et al, 2012), Burkholderia 
might have a pivotal role in the degradation of various insec- 
ticides. To clarify this point, it would be of great interest to 
also employ a deep sequencing survey of 16S rRNA genes 
for various types of organophosphorus insecticide-contaminated 
soils. 



OTHER BACTERIAL GROUPS RESPONDING TO MEP-SPRAYING 

Previous studies have isolated and characterized MEP-degrading 
bacteria from genera Cupriavidus, Pseudomonas, Sphingomonas, 
Corynehacterium, Arthrobacter, and Burkholderia (Tago et al., 
2006; Zhang et al, 2006; Kim et al, 2009). In addition to 
Burkholderia, members of genera Dyella, Ralstonia, Pandoraea, 
and Achromobacter were also isolated as MEP-degrading bac- 
teria in our culture-dependent experiments (Figure 2). Hence, 
this is the first record of MEP-degrading strains from these four 
genera, suggesting that broader taxonomic groups could con- 
tribute to MEP-degradation in natural soil environments. Some 



Frontiers in Microbiology | Systems Microbiology 



August 2014 | Volume 5 | Article 457 | 8 



Itoh et al 



Microbial succession under insecticide-spraying 



86 



Burkholderia sp. FDS-1 AY550913 ^ 
Burkholderia sp. SFA4 AB232336 
Burkholderia sp. MDT2 AB665358 
OTU_N3p2_02 
OTU_N3p1_02 
Burkholderia sp. KM-A AB622644 
Burkholderia sp. KM-G AB622646 
Burkholderia sp. IFA2 AB232326 -jf 
Burkholderia sp. MDT50 AB665361 
★1 



OTU N3 isolates 01 



OTU_N3p2_01 
OTU_N3p1_01 
OTU_S3p2_03 
OTU_S3p1_02 



SBE 
clade 



Dominant OTUs in 
MEP-sprayed soil N 

Burkholderia sp. KFA2 AB232331 

Burkholderia sp. RPE64 AB558208 

Burkholderia sp. Y123 CP003087 

Burkholderia sp. NF100 AB025790-^- 

Burkholderia sp. TFA1 AB232337 -Jf 

Burkholderia sp. TFA2 AB232338 

Burkholderia sp. SFA1 AB232333 

' B. andropogonis X67037 



:fi! 



S. glatheiY 17052 



B. multivorans Y18703 
S. gladioli NR 102847 

S. vietnamiensis AF097534 
S. fetens AM747628 



OTU. 


_S3p1_01 




OTU 


_S3p2_01 




OTU 


S3 isolates 


01 


OTU. 


_N3p1_03 




OTU 


N3 isolates 


02 



> 



Dominant OTUs in 
MEP-sprayed soil S 



B. cepacia AF31 1972 
B. glumae DQ355167 

OTU_S3p1_03 

B. pseudomallei U91 839 



Bcc 
clade 



83 



OTU S3 isolates 03 



B. metallica AM747632 
Burkholderia sp. KFA1 AB232330 
Burkholderia sp. SFA3 AB232335 
Burkholderia sp. SFA2 AB232334 -jc 
B. caribensis AY586519 

OTU_N3p2_03 
OTU_S3p2_02 

S. sordidicola AF5 12827 



PBE 

clade 



71 B. phenazinium U96936 
B. fungorum AF215705 
6S| I a phytofirmans AB741457 

B. kururiensis AB024310 

B. unamae AY221 956 



Pandoraea pulmonicola AF1 391 75 



FIGURE 6 | Phylogenetic analysis of Burkholderia strains dominating 
in MEP-sprayed soils. A neighbor-joining tree inferred from aligned 
255 bp sequences of 16S rRNA gene is shown. Phylogenetic placements 
of OTUs (OTUs at 100% identity thresholds) of Burkholderia, identified 
from the MEP-degrading isolates and deep sequences, were estimated. 
Among dominant OTUs in each soil samples, top 3 OTUs (see Figure 5) 
were included in the analysis. Red and blue colored OTUs are detected 
from soil S and soil N, respectively. Outline letters on a filled background 
indicate OTUs determined by the isolation. Colored letters on a white 



strains oiRalstonia and Achromobacter have already been reported 
to degrade organophosphorus compounds that are analogous 
to MEP, particularly fenamifos (organophosphorus nemati- 
cide) and methyl parathion (organophosphorus insecticide), 
respectively, (Zhang et al., 2005; Cabrera et al., 2010). Similar to 



background indicate OTUs detected by the deep sequencing. Asterisks 
indicate MEP-degrading Burkholderia strains. Clade names of SBE 
(stinkbug-associated beneficial and environmental group), Bcc 
{Burkholderia cepacia complex), and PBE (plant-associated beneficial and 
environmental group), respectively, according to Itoh et al. (2014), Coenye 
et al. (2001), and Suarez-Moreno et al. (2012), are shown on the right 
side. Bootstrap values higher than 50% are depicted at the nodes. The 
16S rRNA sequence of Pandoraea pulmonicola (AF139175) was used as 
an outgroup. 



Burkholderia strains, these bacteria probably have an ability for 
cross-acclimation toward other organophosphorus compounds. 

In addition to the complete degradation of MEP by a 
single bacterium, the cooperative degradation of MEP has 
also been reported in vitro. In this case, MEP-hydrolyzing 
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and MEP-hydrolysate-degrading bacteria cooperate for complete 
consumption of the chemical compound (Katsuyama et al., 
2009). Previous studies have shown that some Burkholderia 
strains, though ineffective against MEP, are able to degrade 
its hydrolysate, 3-methyl-4-nitrophenol (Kim et al, 2007a; 
Katsuyama et al., 2009). Hence, some of the increasing bac- 
teria identified by deep sequencing (Figure 4) from the MEP- 
sprayed soils might be directly involved in the hydrolysis of MEP, 
while others might contribute to the degradation of the result- 
ing hydrolysates. Such a synergistic interaction among bacterial 
species could promote efficient degradation of MEP in soils. 

METHYL0TR0PHS COULD PLAY A ROLE IN MEP-DEGRADATION 

Notably, deep sequencing revealed that the frequency of 
Methylobacillus, Methylobacterium, and Methylophilus, which 
were never detected by the culture-dependent approach in this 
study and from previous researches (Tago et al., 2006; Zhang 
et al., 2006; Kim et al, 2009), also dramatically increased in 
soils after MEP-spraying (Figure 4). Most members of these gen- 
era are methylotrophic, that is, they actively use one-carbon 
compounds such as methanol as a sole carbon source (Kolb 
and Stacheter, 2013). Past studies have demonstrated that bacte- 
ria break down MEP into 3-methyl-4-nitrophenol and dimethyl 
phosphate during the initial hydrolysis step (Spillner et al., 1979; 
IPCS, 1992) (Figure 7). Though the succeeding stages for biode- 
grading the former aromatic compound have been described in 
detail (Spillner et al, 1979; IPCS, 1992; Hayatsu et al, 2000), 
the fate of dimethyl phosphate in soil still remains unclear. In 
the case of the degradation of paraoxon, an analog of MEP, 
the initial hydrolysis process produces diethyl phosphate, which 
is then further hydrolyzed into two molecules of ethanol by 
Delftia acidovorans and Pseudomonas putida through phosphodi- 
esterase and phosphatase enzymes (Tehara and Keasling, 2003). 



Although speculative, dimethyl phosphate may also be hydrolyzed 
into two molecules of methanol by a similar route during MEP- 
degradation. If so, the methylotrophs detected in this study 
may consume and utilize such generated methanol for growth 
(Figure 7). 

It has been reported that accumulated methanol is harmful to 
organisms that are unable to metabolize it (Martin-Amat et al., 
1978). In fact, growth of MEP-degrading Burkholderia strains iso- 
lated in this study was depressed by a tiny amount of methanol 
contamination (Figure S4). If methylotrophs rapidly consume 
any methanol generated during MEP-degradation, they might 
protect other bacteria involved in the overall degradation pro- 
cess from the toxicity of this by-product. This could probably 
lead to better growth of such relevant bacterial communities, 
thus contributing to more efficient and complete mineraliza- 
tion of MEP. This hypothesis could be verified through future 
evaluations on the (1) degradability of dimethyl phosphate by 
methylotrophs and (2) co-culture of MEP-degrading bacteria and 
methylotrophs. 

CONCLUDING REMARKS 

This study demonstrated the effect of MEP-spraying history 
to the microbial succession under MEP application and iden- 
tified key bacteria for MEP-degradation by performing both 
culture experiments and high-throughput sequencing. Owing to 
the enormous information generated by the deep sequencing of 
bacterial 16S rRNA gene, we comprehensively revealed the com- 
munity structures of bacterial species responding to MEP, and 
found previously unseen members, methylotrophs, that may play 
an important role in the complete degradation of MEP in soil 
environments. In contrast, it should be noted that the partial 
sequencing of 16S rRNA gene gave us only indirect informa- 
tion about the MEP-degrading strains. In fact, our previous study 
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FIGURE 7 | Proposed biodegradation pathway of MEP In addition to the 
utilization of 3-methyl-4-nitrophenole by Burkholderia and other bacteria, 
methanol generated from hydrolysis of dimethyl phosphate might be used as 
a carbon source by soil bacteria like methylotrophs. + P=S was rapidly 



substituted to P=0 by oxidation in environment (Bajgar, 2004; Karpouzas and 
Singh, 2006). +t(CH 3 0) 2 POOH + 2H 2 0 -> 2CH 3 OH + PO^ + 3H+. ttt The 
ring cleavage process has been reported in the biodegradation of MEP or 
MEP analog (Hayatsu et al., 2000; Zhang et al., 2006). 
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reported that even Burkholderia strains exhibiting high (>99.8%) 
sequence similarities in 16S rRNA gene possess quite differ- 
ent MEP-degrading activities (Kikuchi et al., 2012). Sequencing 
of functional genes involved in the MEP-degradation pathway 
might give more direct information for MEP-degradation in 
environmental soils, although MEP-degradation genes are not 
well understood (Singh, 2009). In addition, fungal biodegrada- 
tion of MEP has been also reported (Baarschers and Heitland, 
1986), although here we focused only on bacteria. Metagenomic 
approaches of functional genes and more comprehensive sur- 
vey including fungal communities in soils should improve our 
knowledge of MEP-degradation by soil microorganisms. 
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